Squeezed single-atom laser in a photonic crystal 
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We study non-classical and spectral properties of a strongly driven single-atom laser engineered within a 
photonic crystal that facilitates a frequency-dependent reservoir. In these studies, we apply a dressed atom 
model approach to derive the master equation of the system and study the properties of the dressed laser under 
the frequency dependent transition rates. By going beyond the secular approximation in the dressed-atom cavity- 
field interaction, we find that if, in addition, the non-secular terms are included into the dynamics of the system, 
then non-linear processes can occur that lead to interesting new aspects of cavity field behavior. We calculate 
variances of the quadrature phase amplitudes and the incoherent part of the spectrum of the cavity field and show 
that they differ qualitatively from those observed under the secular approximation. In particular, it is found that 
the non-linear processes lead to squeezing of the fluctuations of the cavity field below the quantum shot noise 
limit. The squeezing depends on the relative population of the dressed states of the system and is found only if 
there is no population inversion between the dressed states. Furthermore, we find a linewidth narrowing below 
the quantum limit in the spectrum of the cavity field that is achieved only when the secular approximation is 
not made. An interpretation of the linewidth narrowing is provided in terms of two phase dependent noise 
(squeezing) spectra that make up the incoherent spectrum. We establish that the linewidth narrowing is due 
squeezing of the fluctuations in one quadrature phase components of the cavity field. 

PACS numbers: 42.55.Tv, 42.50.Ct, 42.50.Ar 



I. INTRODUCTION 



A single-atom laser, in which no more than one atom is 
present in an optical resonator, is one of the key systems to 
study quantum effects in the interaction of electromagnetic 
fields with matter lH 12. S B S • Although a single-atom laser 
is admittedly an elementary model, it has the advantage over 
a multi-atom laser that in any practical realization of the laser 
one is not concerned with many difficulties such as fluctua- 
tions of the number of atoms. With the present trapping and 
cooling techniques, the atom can easily localized in a small 
region within the Lamb-Dicke regime. Because of its simplic- 
ity, many quantum features of single-atom lasers, such as sub- 
Poissonian photon statistics, photon antibunching, squeezing 
and vacuum Rabi splitting have been predicted and experi- 
mentally observed 161 17|]. 

The primary obstacle in the realization of one-atom lasers is 
spontaneous emission. The reason is that spontaneous emis- 
sion is a source of noise that leads to emission of photons 
into modes different from the cavity mode. Therefore, it is 
not surprised that many schemes have been proposed to re- 
duce spontaneous emission. It has been demonstrated that 
spontaneous emission can be reduced if an atom is located 
in a squeezed vacuum |8], or inside a cavity that facilitates 
an exclusive spatial selection of radiation modes coupled to 
the atom 191 UOl [ill [l2ll . The cavity appears as a frequency 
filter for spontaneous emission. A further reduction of spon- 
taneous emission can be achieved through a dynamical means 
that unlike the conventional method of coupling an atom into 
the cavity mode, one could first drive the atom with a strong 
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laser field and then cotiple the resulting dressed-atom system 
with the cavity mode 11131 \\M . When the Rabi frequency of 
the driving laser is much larger than the cavity bandwidth, 
one can tune the cavity mode to one of the dressed-atom tran- 
sition frequencies thereby eliminating spontaneous emission 
on the other transitions. Within this approach, a wide variety 
of quantum and spectral features, such as atomic population 
inversion |15, 16, 17], and dynamical suppression and nar- 
rowing of the spectral lines lllSl [1911 has been predicted and 
some of them verified experimentally |[20l[2ll[22l[23[l . 

Recently, Florescu et at. uM have appHed the dressed-atom 
approach to a single-atom laser that incorporates the dynam- 
ical suppression of spontaneous emission inside a photonic 
crystal. They have applied the secular approximation in the 
coupling of a dressed two-level atom to the cavity mode that 
the Rabi frequency of the dressing field is much larger than 
the coupling constant between the atom and the cavity field. 
By appropriate tuning of the cavity mode and the dressed- 
atom transition frequencies, they have shown that the vari- 
ances of the cavity field amplitude and the linewidth of the 
cavity field spectrum can be reduced to the quantum shot noise 
limit. These effects result from the filtering property of the 
photonic crystal that effectively forms a frequency dependent 
reservoir of a specific spectral function (step function) of the 
radiation modes. A photonic crystal that is a periodic dielec- 
tric structure, can prohibit light propagation over a continuous 
range of freq uencies, irrespective of the direction of propaga- 
tion 1251 [2611 . It is also known that in photonic crystals ex- 
tremely small micro-cavity mode volumes and very high cav- 
ity Q factors can be realized 127, 28, 291 . It is associated with 
the unique properties of photonic crystals, i.e., the photonic 
density of states within or near a photonic band gap can al- 
most vanish or exhibit discontinuous changes as a function of 
frequency with appropriate engineering, which is essentially 
different from its free-space counterpart. 



Although the secular approximation is often entirely ad- 
equate to describe the interaction of the cavity field with a 
dressed-atom system, there are circumstances when the non- 
secular terms might be important, for example in the strong 
coupling limit of the cavity field to the dressed-atom system. 
One could argue that even under these circumstances, the non- 
secular terms would only make small corrections to the dy- 
namics of the system. In the following, we show that this is 
not the case, they may lead to new aspects of cavity field be- 
havior that differ qualitatively from those observed under the 
secular approximation. To show this, the dressed-atom ap- 
proach for a driven two-level atom placed inside a photonic 
crystal is generalized to deal with non-secular terms in the 
atom-cavity field interaction. The method used in this paper 
is similar to that used by Florescu et al. [24] but with one sig- 
nificant difference, we do not make the secular approximation 
in the couphng of the cavity mode to the dressed-atom system. 
Instead, we consider the role of the non-secular terms in the 
operation of a strongly driven single-atom laser We present 
solutions for the steady-state variance of the quadrature com- 
ponents of the cavity field amplitude and the incoherent part 
of the spectrum of the cavity field. In the course of the calcu- 
lations, we observe that the inclusion of the non-secular terms 
results in a double coupling of the cavity field to the dressed 
atom and leads to non-linear terms in the master equation of 
the system that reflect a possibility of non-linear processes be- 
tween the dressed states when the secular approximation is not 
made. Our results demonstrate that these non-hnear processes 
lead to the reduction of the fluctuations of the cavity field am- 
plitude below the quantum shot noise limit. In addition, we 
find that the spectral line of the cavity field may now be nar- 
rowed below the quantum limit giving a subnatural linewidth 
of the emitted field. A qualitative understanding of the ori- 
gin of the line narrowing is obtained by writing the incoherent 
part of the cavity field spectrum in terms of the phase depen- 
dent noise (squeezing) spectra. We show that the narrowing 
results from negative values of the squeezing spectrum of one 
quadrature phase component of the field. 

The paper is organized as follows. In Sec. II, we present 
the model Hamiltonian and derive the master equation for a 
strongly driven two-level atom coupled to a single-mode cav- 
ity and placed inside a photonic band gap material that acts 
as a frequency dependent reservoir. In Sec. Ill, we derive 
the basic equations of motion for the expectation values of the 
atomic and field correlation functions are find their steady- 
state values. We then apply them in Sec. IV to investigate 
variances of the quadrature phase components of the cavity 
field that determine the fluctuations of the cavity field ampli- 
tude, and in Sec. V to the calculations of the spectrum of the 
cavity field. We calculate the incoherent part of the spectrum 
and then analyze its properties in terms of the squeezing spec- 
tra. Finally, we summarize the results in Sec. VI. 



11. MASTER EQUATION 

We consider a single two-level atom with ground state |1) 
and excited state |2} separated by the transition frequency tOa- 



The atom is coupled to a single mode of a high-Q microcavity 
engineered within a photonic crystal with coupling constant g, 
and is driven by a coherent external laser field of a frequency 
lol and the resonant Rabi frequency e. In addition, the atom is 
damped at the rate 7 by spontaneous emission to modes other 
than the cavity mode. In practice this model can be realized by 
embedding a quantum dot in a dielectric microcavity (defect) 
placed within a two-mode waveguide channel in a 2D PBG 
microchip [24]. In a photonic band gap material, one mode of 
the waveguide channel is engineered to produce a large dis- 
continuity in the local photon density of states near the atom, 
and another mode is used to propagate the pump beam. By 
suitable engineering, it is possible to realize a strong coupling 
of the quantum dot to both the pumping waveguide mode and 
the high-Q cavity mode [24]. For simplicity, we treat the driv- 
ing external field classically and work in the interaction pic- 
ture. The total system is described by a Hamiltonian that un- 
der the electric-diople and the rotating-wave approximations 
can be written as 



H — Hf) + Hi, 



(1) 



where the first term 



Hq = hAca^a + -hAaas 



+he{ai2 + a2i) + h^ Axa\ax (2) 



is the noninteracting Hamiltonian of the driven atom plus 
the cavity mode plus the photonic crystal radiation reservoir 
modes, and the second term 



Hi = ihg{a' ai2 - (J2ia) 

+ih'^gx{uj\){a\ai2 ~ (J2ia\) 

X 



(3) 



is the interaction Hamiltonian between the atom and the cavity 
mode and the photonic crystal vacuum radiation modes. 
Here, a and at are the cavity-mode annihilation and cre- 
ation operators, a\ and a^ are the photonic crystal radia- 
tion reservoir annihilation and creation operators, crij are the 
bare atomic operators, aij — \i){j\ {i,j — 1,2), and 0-3 = 
f22 — cTii describes the bare atomic inversion. The parameter 
Aa — LUa — lul dcnotes the detuning of the atomic resonance 
frequency uta from the frequency cj^ of the driving laser field, 
Ac = ojc — ^L is the detuning of the cavity frequency from 
the frequency of the laser field, and A\ = uj\ — ujl is the 
detuning of the reservoir frequency uj\ of a mode A from the 
laser frequency. 

The coefficient g describes the strength of the coupling be- 
tween the atom and the cavity mode that we assume to be 
constant independent of frequency, and g\{LU\) describes the 
strength of the coupling between the atom and the vacuum 
modes of the photonic crystal. It contains the information 
about the frequency dependent mode structure of the photonic 
crystal and can be written as 



9\i(^x) ^ 9\D{ujx), 



(4) 



where g\ is a constant proportional to the dipole moment of 
the atom, and D{u!x) is the transfer function of the photonic 
crystal, the absolute value square of which can be identified 
as the Airy function of a frequency dependent radiation reser- 
voir L30J. 

The explicit form of D{u!\) depends on the type of ra- 
diation reservoir used. For a photonic band gap material, 
the transfer function is in the form of the unit step function 
\D{u}\)\^ — u{uj\ — LOb), where Wf, is the photonic density of 
states band edge frequency. Thus, |£'(wa)P = for lox < ujb 
and |£'(wa)P = 1 for ijJ\ > cjf,. Since the frequencies 
uJx < ^b are forbidden in the band gap material, it is pos- 
sible to selectively eliminate spontaneous emission at some 
frequencies of the driven atom. 

The strong driving field can be viewed as a dressing field for 
the atom. Therefore, we begin by diagonalizing the atomic 
part of the Hamiltonian together with the interaction of the 
atom with the laser field 



Haf = -hAaaa + he{ai2 + Cr2l), 



(5) 



to find the eigenstates (dressed states) of the combined atom 
plus driving field system. Since the driving field is treated 
classically in our calculations, we find the so-called semiclas- 
sical dressed states 



|i) = cos(/)|l) +sin</)|2), 
|2) = sin0|l} -cos</>|2), 



(6) 



where cos^ (p — {I + Aa/57)/2 with the angle defined such 
that < (f> < tt/2. The dressed states form non-degenerate 
doublets that are separated in energy by huj^, and the states of 
the doublet ai-e split by h2Q, where O = (4e^ + A^y/^ is the 
Rabi frequency of the detuned field. 

We now couple the dressed states to the cavity field and to the 
photonic crystal vacuum modes. First, we replace the atomic 
operators by the dressed-state operators 



cri2 = -- sin(20)i?3 + sin^ 0i?2i - cos^ (t)Ri2, 

o"2i = — - sin(20)i?3 + sin^ 0i?i2 — cos^ (f>R2i, 
as = -cos(2(/.)i?3+sin(2(/))(i?i2 + i?2i), (7) 

where Rij = \i) {j \ are the dressed-atom dipole operators and 
^3 = R22 — Rii- Next, we perform the unitary "dressing" 
transformation of the interaction Hamiltonian 



with 



Hi = eyip{iHot)Hi exp(— liJo^), 



Ho = flRs + Aca'^a + ^ Axa{ax, 
A 



(8) 



(9) 



and obtain the interaction Hamiltonian between the dressed 
atom and both the cavity mode and the photonic crystal vac- 



-2n)t 



uum modes 

Hi = i% rscati^ge'^-* + c2ati?i2e'^'^'=" 
-s2ati?2ie^(^=+'")* - H.c.) 



,^a{R2ie'<^^'+^''^* - H.c.) , 



(10) 



where s = sin and c = cos <p. 

The Hamiltonian (fTol l describes the interaction of the dressed 
atom with the cavity field and with the vacuum modes. We see 
that in the dressed-atom picture, the cavity frequency and the 
vacuum modes are tuned to the dressed-state transitions that 
occur at three characteristic frequencies, Ac and Ac ± 251. By 
matching the cavity field frequency to one of the dressed states 
frequencies, we may manipulate the strange of the interaction 
between the driven system and the cavity mode. 

Our aim is to derive from the Hamiltonian ( fTol i the master 
equation for a reduced density operator of the driven atom and 
the cavity field. It is obtained by tracing the density operator 
of the total system over the photonic crystal radiation reservoir 
variables. On carrying out this procedure, it is found that in 
the dissipative part of the master equation certain terms are 
slowly varying in time while the others are oscillating with 
frequencies 2i7 and 451. Since we are interested in the case 
where the Rabi frequency ft is much larger than the atomic 
and cavity damping rates 



rj >7,K, 



(11) 



we can invoke the secular approximation that consists of 
dropping these rapidly oscillating terms. These terms, if kept 
in the master equation, would make corrections to the dynam- 
ics of the system of the order of 7/il, and thus completely 
negligible. We therefore find that after discarding the rapidly 
oscillating terms in the dissipative part of the master equation, 
the time evolution of the reduced density operator is of the 
form 



^ = gsc[a^ R^e'^^' - aR. 
dt ^ ^ 



:3e ^^=Sp] 



+ 5c2[ati?i2e'(^=-20)t _ aR2ie-'^'^^-^'''^\ p] 

- gs2[ati?2ie'(^=+2")* - ai?i2e-^(^=+2")*,p] 

+ LaP + C,p, (12) 



where 



^aP = IpQ {R3PR3 - P) + xT- {R21PR12 - R12R21P) 

+ 27+ {R12PR21 - R21R12P) + H.c. 
CcP ~ —K (2apa' — a^ap — pa^aj , 

describe spontaneous dynamics between the dressed states of 
the system and of the cavity mode, respectively. The parame- 



ters 



70 = s^c^-t\D[uL)? + [^-s^)-ij,, 
7_ = sSl^K-2r!)|2 + 4s2c27p, 
7+ = cS|^(^L + 2r!)|2 + 4sV7p 



(13) 



determine the damping rates between the dressed states of the 
system. They also include the contribution of the dephasing 
rate 7p which may arise from scattering of phonons of the host 
crystal on the atom embedded in the solid part of the dielec- 
tric material. The coefficient 70 corresponds to spontaneous 
emission occurring at two transitions of the dressed atom; One 
from the lower dressed state |1) to the lower dressed state of 
the manifold below and the other from the upper dressed state 
1 2) to the upper dressed state of the manifold below. These 
transitions occur at frequency lol- The coefficient 7+ corre- 
sponds to spontaneous emission from the upper dressed state 
to the lower dressed state of the manifold below and occurs at 
frequency ujl + 217, whereas the coefficient 7_ corresponds 
to spontaneous emission from the lower dressed state to the 
upper dressed state of the manifold below and occurs at fre- 
quency lol — 2ri. 

An essential feature of the master equation (fT2l i is the pres- 
ence of the oscillatory terms in the coherent part of the evolu- 
tion involving the dressed atom and the cavity field. In exper- 
imental practice it might imply that there is no restriction on 
the relation between the Rabi frequency fl and the coupling 
constant g, and the master equation can be applied to both 
situations, where the coupUng constant g is much smaller or 
comparable to the Rabi frequency. 

As we have already mentioned, we work in the strong- 
coupling regime of fi ^ g ^ 7, k. The latter inequality 
has been invoked into Eq. (fT2l l to neglect the effect of non- 
secular terms in the dissipative part of the master equation, 
but to keep the non-secular terms in the interaction of the cav- 
ity field with the dressed atom. This condition is satisfied 
in a typical band-gap material, where strong couplings g and 
significant reductions of the spontaneous emission rate 7 are 
easily achieved due to the confinement of the atom to an ex- 
tremely small volume with a significantly suppressed density 
of the vacuum modes. For example, given the electric dipole 
moment of an atom of // « lO^^^C m^^ and the cavity mode 
volume of y « 10~^ m'^, then the atom cavity coupling will 
be g « 10^'' Hz ll24ll27ll . With the spontaneous emission rate 
of the atom in a free space of 7 « 10* Hz, the inequality 
g ^ 7 is easy satisfied. The parameter k, that describes the 
damping of the cavity field, is equal to zero in an ideal cavity. 
However, in a realistic structure k 7^ 0, that arises from a weak 
coupling of the cavity mode to the waveguide mode of the 
band-gap material. This coupling may result from defects or 
disorders in the waveguide channel caused by the manufactur- 
ing process that can be kept at a minimal level 11241 l2a, [27112911 . 
Effectively, the waveguide channel acts as an external reser- 
voir of vacuum modes to which photons can escape from the 
cavity mode. Thus, the inequality g ^ 7, k is realistic and 
should be easy to satisfy within the current band-gap material 
technology. 

In the following, we focus on the case of Ac « 0, i.e., the 



cavity field tuned close to the central frequency of the dressed 
atom. In this limit, the master equation ( fT2] ) contains terms 
that are time independent and thus corresponding to the res- 
onant interaction of the cavity field with the dressed atom. It 
also contains terms that have an explicit time dependence of 
the form exp(±2if7f). These terms correspond to a dispersive 
(non-resonant) interaction of the cavity field with the dressed 
atom, which induces interesting new property as we will show 
below. 

We first perform a canonical transformation of the master 
equation ( fT2] l by 



exp(iAca^at)pexp(— iAcfl'at), 



(14) 



and find that the master equation for the transformed density 
operator takes the form 



^ = -^[H^ + H'{t),p] 



Cap + -Cc/5, 



(15) 



where 



Hq = gsc {a* Rz — aR^) + Aco'a, 
H'it) = ig [{c^a'^ + s^a) R^e' 



i2nt 



H.c. 



(16) 



Note that the Hamiltonian H'{t) disappears under the secular 
approximation on the interaction of the dressed-atom with the 
cavity field. Due to the presence of the rapidly oscillating 
terms, the Hamiltonian H'{t) can be treated as a perturber to 
the time-independent Hamiltonian Hq. Since g ^ fl, we can 
perform a second-order perturbation calculations with respect 
to g and find an effective Hamiltonian 13 ill that can be written 
as 



^cff = -* 






c") {R3a^a+-R3 



In this case, the master equation (fT2b takes the form 



(17) 



dp 

dt 



9I 



a' a, p 



.91 



4[^3,P] 



CaP + CcP, 



where 



91 = 9SC 



and 



92 



= .9V(^ 



(18) 



(19) 



are the "effective" coupling constants of the cavity field to the 
dressed-atom system. 

The important new feature of the master equation ( fTSl l is in 
the presence of an additional non-linear terms that results in 
an effective double coupling of the cavity field to the dressed 
atom. The first term on the right-hand side of Eq. (fTSl l rep- 
resents a shift of the cavity frequency. It includes a dressed- 
atom inversion dependent shift that arises solely from the non- 
secular terms in the dressed-atom picture. The second term 



represents a shift of the dressed-atom frequencies that also 
arises from the presence of the non-secular terms. The third 
term represents a resonant coupling of the cavity field to the 
dressed atom at frequency lOc ~ lol, and finally the fourth 
term describes a dispersive coupling through the off-resonant 
Rabi sideband frequencies and is determined by non-linear 
two-photon absorption and emission processes between the 
dressed states. Thus, the presence of the non-secular terms 
results in the non-linear coupling of the cavity field to the 
dressed atom and shifts of the cavity and dressed-atom res- 
onance frequencies. The amplitude of these contributions is 
of order g^ /Vl and therefore is small, but not necessarily too 
small to make detectable contributions to the dynamics of the 
system. In what follows, we shall show that the non-linear 
terms in the master equation dTsl l play in fact an important 
role in the dynamics of the system and will study in details 
non-classical and spectral properties of the out-put field of the 
laser operating with a resonator engineered in a photonic band 
gap material. 



III. EQUATIONS OF MOTION 



In this section we investigate the influence of the non-secular 
terms on the properties of a driven single-atom laser The 
properties can all be expressed in terms of single-time and 
two-time expectation values of the dressed-atom and the cav- 
ity field operators. The master equation ( fTSl l enables us to 
derive equations of motion for expectation values of an arbi- 
trary combination of the atomic and cavity field operators. In 
particular, for the dressed-atom population inversion and the 
cavity field amplitudes, we find the following closed set of 
equations of motions 



ter equation leads to the following set of equations of motion 

d 



dt 



{a^a) = - K{a^a) + gi {{R^a^) + {Riu)) 
- i^ {{Rsa^') - {Rsa')) , 



-^(i?3at2) ^ j|i _ (^, + « + 2zA,) (i?3«^') + 2.gi(a^) 
dt \l 

-72(«^^)+^ff(at2)+2^|i(ata), 



dt 



{RWo) = - (71 + K)(i?3a^a) + gi{{a)) + (a)) 

-I2{a^a)-i^l{{a'^^)-{a')), (21) 



and equations of motion for {a?) and {R-^a?) are obtained by 
the Hermitian conjugate of the equations of motion for (a^^) 
and (_R3a^^), respectively. 

The quantities 71 and 72 occurring in Eqs ( l20l l and ( I2TI) are 
given in terms of the damping rates 7+ and 7_ , and are defined 
as 



2„2 



71 = 7+c +7_s + 87pS c 

72 = 7+c^ -7_s''. 



(22) 



We see that the quantities 71 and 72 depend only on the tran- 
sition rates centered at the dressed states resonances corre- 
sponding to the Rabi sidebands. Moreover, the quantity 72 
can be negative, whereas 71 is positive for all values of the 
parameters involved. In addition, 71 > 72 even in the absence 
of the dephasing process, 7p = 0. 

IV. STEADY-STATE SOLUTIONS 



— {i?3} =-72 -7i(i?3), 

^(a)=- QK-zA,j(a)+.gi(i?3) 

— (i?3a) =51 - f 71 + -K-iAcj (i?3a) 
^ ' 2Vt^ ' D.^ ' 



(20) 



and equations of motion for (a^) and (Rsa^) are obtained by 
Hermitian conjugate of the above equations. 

For the expectation values involving higher order combina- 
tions of the operators, such us the number of photons, the mas- 



We proceed here to discuss the steady-state properties of the 
correlation functions for Ac — that are listed in Appendix 
I. In this simplified case, we obtain solutions in a physically 
transparent form that will allow one to gain physical insight 
into how the non-secular terms and structured band gap mate- 
rial can modify the properties of the dressed-atom system and 
the cavity field. The exact steady-state solutions of Eqs. ( |20l ) 
and jTH . valid for an arbitrary detuning Ac are quite lengthy 
and will be computed numerically in the next section. 

Before commencing our analysis of the properties of the cav- 
ity field, we first look into properties of the dressed atom. We 
observe that the equation of motion for the expectation value 
of the atomic operator R3 is decoupled from the remaining 
equations. Hence, it has a simple steady-state solution 



{R3) = -^ 



72 
71' 



(23) 



This shows that the steady-state value of the dressed inversion 
operator is not affected by the non-secular terms. However, 



it depends crucially on the quantities 71 and 72 that, on the 
other hand, depend on the mode structure of the band gap ma- 
terial. The inversion varies between + 1 and — 1 corresponding 
to locking the population in the dressed state |2) and |1), re- 
spectively. According to Eq. ( l22b . it happens when there is 
no spontaneous emission on one of the Rabi frequency of the 
dressed system. 

Let us now examine the stationary properties of the cavity 
field. In the steady-state regime and for the case of the cavity 
frequency and the atomic transition frequency tuned on reso- 
nance with the central frequency of the dressed-atom system, 
Ac — and A^ — 0, the expectation values of the cavity field 
operators can be expressed as 



7?(7?-7l) 




K725 2 


2(7?-7l)l 


71 (271 + k) 




1- 


(7?-7l) 



K72 (271 + k)Q 



(3«; + 47i)(72-72) 



(k + 7i)(k + 271)2 



71 (71 + k) 



" (7i - 7I) 


\ 


71 (71 + 1^) (271 + '«) 


/• 


2(7?-7l) , gf [^ 
71 (271 + «) 2172 


71 (71 + k) 


u (7? - 7I) 


}), <2„ 


71 (71 + f^f (271 + f^f 



where m = (71 + k) (271 + k) (271 + 3k) + (471 + 3k) 72. 
Similarly to the dressed-atom inversion, the properties of the 
expectation values of the cavity field depend on the parameters 
71 and 72 and thus are also sensitive to the mode structure of 
the band gap material. A particular situation of interest is that 
of 71 = 72 which, according to Eq. (|22]) occurs when the 
dephasing rate jp — and the lower Rabi sideband frequency 
is adjusted to be inside the forbidden frequency region of the 
band gap material, 7_ = 0. As we see from Eq. dZST l. the 
cavity field correlation functions are then enhanced due to the 
combined effect of the band gap material and the non-secular 
processes. 

The non-secular terms not only enhance the number of pho- 
tons, but also lead to a shift of the cavity resonance frequency. 
This is shown in Fig.lT] where we plot the steady-state expec- 
tation value of the number of photons in the cavity field as a 
function of the detuning Ac in the band gap material configu- 
ration 7_ = 0, in which spontaneous emission is forbidden at 
the lower Rabi sideband frequency. The graphs are computed 
from the exact steady-state solution of Eq. ( I2TI 1 that is pre- 
sented in the Appendix I. It is evident that in the presence of 
the non-secular terms the number of photons attains the maxi- 
mum value in the vicinity of the detuning Ac = — O.I67. The 
shift of the maximum towards the negative detuning Ac indi- 



cates a pulling of the oscillation frequency of the cavity field 
away from the passive cavity frequency uJc towards the lower 
Rabi sideband frequency that is inside the forbidden frequency 
region of the band gap material. This is another noteworthy 
feature of the non-secular terms in the dressed-atom picture. 

We note here that the values of the parameters used to plot 
Fig.[T]are consistent with the good cavity limit assumed in the 
derivation of the master equation ( fTSl l and are experimentally 
realistic. We have discussed the validity of the good cavity 
limit in Sec.|II]for realistic parameters of our system, and have 
shown that a small cavity damping and a large atom cavity 
coupling g can be readil y ac hieved within the current band- 
gap material technology ll24ll27ll . 



V. FLUCTUATIONS OF THE CAVITY FIELD 

We now proceed to investigate the fluctuations of the cav- 
ity field amplitude by analyzing the variances of the conju- 
gate quadrature phase amplitudes X^ = ae'^ + a^e^'^ and 
X_ = —i{ae'^^ — a^e^'-^), where 9 is the quadrature phase. 
Previous calculations of Florescu et al. Ii24i1 have shown that 
under the secular approximation and by appropriate tuning of 
the dressed-atom transition frequencies to the band gap ma- 
terial mode transfer structure, the fluctuations of the cavity 
field can be reduced down up to the quantum shot noise limit. 
We now consider the possibility of using the non-secular pro- 
cesses to reduce the quantum fluctuations below the quantum 
limit. There is a good reason to believe that the non-secular 
processes can modify the fluctuations of the cavity field am- 
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FIG. 1 : The steady-state expectation value of the number of photons 
in the cavity field plotted as a function of the detuning Ac in the 
presence (solid line) and in the absence (dashed line) of the non- 
secular terms for il — IOO7, 7+ = IO7, k = O.I7+, g — 10k, 
Aa = 0, and 7_ = 7p = 0. 



plitude. They appear as non-linear processes and it is well 
known that these processes are typical sources used to gener- 
ate squeezed light ^ . 

Let us first examine the normally ordered variance of the 
quadrature component X+ for the simplified case of Ac = 
0. This analysis provides simple analytical formulas for the 
variance and illustrates the role of the non-secular processes 
in the behavior of the cavity field. To determine the variance, 
we make use Eq. dZST l and find that it can be expressed as 



{■.{AX+r ■.)s = Si + S2, 



where 



5i- 



85?(7? 



7l) 



k7^(k + 27i) 



(1 + cos 26*) 



(26) 



(27) 



is the contribution of the linear interaction term, and 

2gh2sin2^r , S2gli-ff - jl){3-f, + 2k) 



52=- 



KJifl 



1 



4^^(1-008 261) 



K7i(k + 7i)(k + 27i)2 
(7?-7l) 



K^i}^ { 7i('« + 7i) 

32ff?(7? - 7l) h7i(^ + 7i) + (571 

K72(k + 7i)2(k + 27i)2 



(28) 
4^)7l] 1 



is the contribution from the non-secular terms. 

Clearly, Si > and 5*2 vanishes under the secular approxi- 
mation. Note that the term S2 is negative for the phase angle 
< 9 < tt/2 and 72 positive. Thus, for an appropriate choice 
of the parameters, the fluctuations of the cavity field can be 
squeezed below the quantum shot noise level. However, when 
we examine Eqs. jZTT i and ( |28] ), we find that the Si term ex- 
ceeds the 5*2 term independent of the values of the parameters 
involved. This means that squeezing is possible only when the 
term 5*1 is suppressed. Inspection of Eq. ( |28] ) shows that the 
Si term is suppressed when 71 =72. According to Eq. (|22]) . 
this conditions can be achieved by suppressing spontaneous 
emission on transitions at the lower Rabi frequency ljl — 2U. 
Thus, the immediate consequence of the fact that 71 — 72 
is that the dominating positive term Si vanishes, leaving the 
variance depending solely on the term 5*2 that arises from 
the presence of the non-linear processes between the dressed 
states. Note the dependence of the "squeezer", the S2 term 
of the variance, on the ratio ^2 111 that according to Eq. ( |23] l 
determines the inversion between the dressed states. Appar- 
ently, for phase angles ^ < Q < tx 12, squeezing is obtained 
only when 72 > 0, i.e. when there is no inversion between the 
dressed states and vice versa, for phase angles 7r/2 < 6* < vr, 
squeezing is obtained for a negative 72, i.e. when there is an 
inversion between the dressed states. Moreover, the optimal 
condition for squeezing is to maintain 71 = 72 . These analy- 
sis indicate that a photonic band gap material can be employed 
to generate a squeezed field from a single-atom laser. 

Figure|2]shows the normally ordered variance (: (AX+)^ :)s 
as a function of the quadrature phase angle d. As the phase 
change, the variance oscillates from negative to positive val- 
ues. Negative values mean squeezing of the fluctuations of 



the quadrature phase amplitude below the quantum shot noise 
level. The maximum negative value of the variance, corre- 
sponding to the maximum squeezing occurs at 6* « 0.55. 
However, the maximum value of squeezing obtained for Ac = 
might not be the greatest value possible due to the shift 
of the cavity resonance. Actually, it appears for a detuning 
slightly shifted from the cavity resonance. This is shown in 
Fig. |3] where we plot the variance as a function of the scaled 
Rabi frequency e/A^, for 71 = 72 and different detunings 
Ac. We see that the variance is negative for all values of the 
Rabi frequency, thus indicating that the cavity field can be in a 
squeezed state irrespective of the strength of the laser field. As 
the Rabi frequency increases, the variance decreases and then 
saturates to a constant value. The saturation level decreases 
with increasing negative detuning Ac and attains the maximal 
negative value, corresponding to the maximum squeezing, for 
Ac = -0.267. 

We may conclude this section that the filtering properties 
of the band gap material and the treatment of the interaction 
between the cavity field and the dressed atom beyond the sec- 
ular approximation produces important modifications in the 
variance of the quadrature field amplitudes, greatly reducing 
the fluctuations below the quantum limit. 



VI. SPECTRUM OF THE CAVITY FIELD 

In this section we investigate the spectrum of the cavity field 
and point out one more difference between the secular and 
non-secular single-atom lasers. It has been shown by Florescu 




LOO 



FIG. 2: The steady-state variance {: (AX-f)^ :)s of the quadra- 
ture component of the cavity field as functions of the quadrature 
phase Q j-K in a full photonic band gap and for vanishing detuning 
Ac = 0. All other parameters are same as in Fig. 1. 
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et al. fM\ that in the secular approximation, the spectral line 
of the cavity field can be narrowed to the quantum limit. Here, 
we will illustrate that the linewidth of the spectrum can be 
further reduced that it may go below the quantum limit if the 
non-secular processes are included into the dynamics of the 
system. 

The stedy-state spectrum of the cavity field is defined as 
the Fourier transform of the two-time correlation function 
{a^{t)a)s of the cavity field operators. In the interaction pic- 
ture, the spectrum takes the form 

/•oo 

S{uj) = 2Re dte'('^-'^^)*(at(i)a),, (29) 

where uj is the spectral frequency, the subscript s represents 
the average over the steady-state values of the field operators, 
and the operators without time argument refer to their steady- 
state values. 

We may write the cavity field operators in the form a{t) — 
{a{t)) + 6a{t), where 6a{t) are the fluctuations operators that 
describe fluctuations of the cavity field about its average value. 
In this case, the spectrum can be decomposed into a sum of 
two terms 



S{UJ) = Sel{i0) + Sin{uj), 

where 

Sei{uj) = 2Tr {a^ ) s {a) sS{uj - ujl) 
is the elastic component of the spectrum, and 

/"OO 

5™(w) = 2Re / dte'(— ^)*(at(i), a). 



0.00 



(30) 
(31) 

(32) 




8/|Aa 



FIG. 3: The normally ordered variance (: (AX+)^ :)s as a function 
of the scaled Rabi frequency e/Aa in the band gap material config- 
uration of 7- =0 for e = 0.8, n = IOO7, 7+ = IO7, K = O.I7+, 
g = 10k, 7p = and different detunings Ac: Ac — (dotted line), 
Ac = — O.I67 (dashed line) and Ac = — O.267 (solid line). 



is the inelastic component with 

(at(t),a), = {a\t)a)s - {a^ {t}) s {a) ^ . (33) 

In order to calculate the two-time correlation function appear- 
ing in E g . (|32] |. we make use of the quantum regression the- 
orem II32I I33II from which it is well known that the two-time 
correlation functions obey the same equations of motion for 
t > as the corresponding one-time correlation functions. 
Therefore, we may used Eqs. ( |20] | and ( 1211 1 to compute the 
spectrum. The spectrum depends, of course, on the state of 
the system which we take to be the stationary state. 

We will compute numerically the spectra using the equations 
of motion ( |20l ) and ( |2TI ). However, in order to gain insight into 
the physics involved, we will present an analytical solution for 
the case of Ac = 0. 

For purposes of numerical computation, it is convenient to 
take the Laplace transform of the equations of motion for the 
correlation functions that transforms them into a set of alge- 
braic equations, and evaluate the spectrum from 

S,n{^) = 2Re[(a^(p),a),]|p^,(^_„^), (34) 

where p is a complex Laplace transform parameter. The ex- 
act solution for the Laplace transform of the correlation func- 
tion with Ac = is given in the Appendix II. We use this 
exact solution to evaluate the spectrum in the presence of the 
non-secular processes. We compare the spectra with those ob- 
tained under the secular approximation where the non-secular 
processes are ignored. The results are show in Figs. |4] and |5] 




-2 2 

(o-coO/y 



FIG. 4: The incoherent part of the spectrum of the cavity field as a 
function of the frequency (uj — ujl)/') calculated with (solid line) and 
without (dashed line) the non-secular terms for 7-/7+ = 0.0001 
and Ac — 0. All other parameters are same as in Fig. 1. 

Let us first compare the spectra obtained with the non-secular 
processes with that obtained under the secular approxima- 
tion. Figure |4] shows the incoherent part of the spectrum of 
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the cavity field in the band gap material configuration with 
7_ — 0.000l7_|_. The spectrum consists of a single line cen- 
tered at frequency to ~ ojl whose the profile depends on 
whether the non-secular processes are or are not included. 
The shape of the line is not exactly a Lorentzian when the 
non-secular processes are included, but becomes a Lorentzian 
when these processes are ignored. One may notice a signif- 
icant linewidth narrowing below the natural width when the 
non-secular processes are included. 

Figure |5] illustrates the effect of the cavity detuning on the 
spectrum. We see that the maximum narrowing of the spectral 
line actually occurs not at Ac = but for the detuning Ac = 
—0.267 that, according to Fig. [3] corresponds to the maximum 
squeezing in one of the quadrature phase components of the 
cavity field. This indicates that the spectrum may be directly 
affected by squeezing in the cavity field. 




(co-coO/y 



FIG. 5: The incoherent part of the spectrum of the cavity iield as a 
function of the frequency {lo — '^Lj/'y for different detunings Ac: 
Ac = (dashed line), Ac = -0.16 (dotted line) and Ac = -0.26 
(solid line). All other parameters are same as in Fig. 1. 

In order to gain insight into the source of the narrowing of the 
spectral line, we take the limit of Ac — Aq = 0, 7_ = and 
derive from the general solution, Eq. ( l39l ). a simple analytical 
formula for the spectrum 



Sin{i^) 



l&ngf 



r22 [k2 + 4(CJ - LOlY 



(35) 



This result is in marked contrast to that obtained under the 
secular approximation ll24ll . In the case of vanishing mode 
density on the lower Mollow sideband (7_ = 0, which corre- 
sponds to a full photonic band gap) and no dipolar dephasing 
(7p = 0), the spectrum of the cavity field consists only of 
the elastic component under the secular approximation Il24ll . 
When the non-secular processes are included, the incoherent 



part of the spectrum is present and shows an interesting prop- 
erty that instead of being a simple Lorentzian, it is instead in 
the form of a squared Lorentzian. 
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FIG. 6: The incoherent part (solid line) of the cavity field spectrum 
together with the squeezing spectra X+ (lo) (dotted line) and X- (lo) 
(dashed line) plotted as a function of (lo — ijJL)l'^i for Ac = — O.267 
and Q = 0.8. All other parameters are same as in Fig. 5. 

It is well known that a squared Lorentzian can be decom- 
posed into a difference between two Lorentzians Il34ll35ll36ll . 
The immediate consequence of one of the Lorentzians being 
negative is to produce a spectrum which fell off as lo^"^ in 
the wings, rather than the lj^^, which would result if both 
Lorentzians were positive. It has been shown that the negative 
weight of one of the Lorentzians can be related to squeezing 
in the field. Thus, the naiTowing of the spectral line, seen in 
Figs.|4]and|5]can be related to squeezing produced in the non- 
linear interaction of the cavity field with the dressed atom. 
Physically, this effect can be best understood and explained 
in terms of the so-called squeezing spectrum. The incoher- 
ent spectrum can be decomposed into the sum of two phase 
dependent squeezing spectra 



5„H = i[x+H + x_H], 



(36) 



where 



J — CXj 

X_(w) = / dte'('^-'^^)*(:X_(t),X_ :), (37) 



are the squeezing spectra of the normally ordered quadrature 
phase components of the cavity field amplitude. Negative val- 
ues in the squeezing spectra indicate squeezed fluctuations of 
the cavity field. 
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Figure |6] shows the incoherent spectrum together with the 
squeezing spectra calculated for the same parameters as in 
Fig. m but with Ac = — O.267 corresponding to the max- 
imum squeezing in the X^ quadrature components of the 
cavity field. We see that the squeezing spectrum X+{uj) is 
negative over the whole range of frequencies. The fact that 
X^{uj) is negative indicates that the X^ quadrature compo- 
nent, that is obtained by an integration of the squeezing spec- 
trum X_|_(w) over the frequency oj, is squeezed. Thus, we 
may conclude that the spectral linewidth narrowing is due to 
the common effect of squeezing and the frequency dependent 
reservoir formed by the photonic band gap material. 



VII. CONCLUSIONS 

We have examined nonclassical and spectral features of a 
single-atom laser without the use of the secular approxima- 
tion for the interaction of the cavity field with a dressed-atom 
system located inside a photonic band gap material. The mate- 
rial appears as a frequency dependent reservoir for the dressed 



atom. The inclusion of the non-secular terms leads to an effec- 
tive double coupling of the cavity field to the dressed states of 
the system. We have found that by a suitable matching of the 
dressed-atom transition frequencies to the guiding frequency 
of the band gap material, it is possible to account for squeez- 
ing and narrowing of the spectral line of the cavity field be- 
low the quantum shot noise limit. These features are achieved 
by a combination of a strongly frequency-dependent reservoir 
and the non-secular processes. The non-linear processes are 
absent under the secular approximation and present a newly 
encountered features of a single-atom laser. 
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APPENDIX I 



and 



In this appendix, we present the exact steady-state solution of 
the equations of motion ( l20l i and (l2Tl i for the case of Ac = 
that includes all non-secular terms up to the order g'^/fl'^. 



{a)s = -2gi[2aFi +72^(271 + K)^]/i^i, (38.2) 

(a^), = {[-2(a - a')F2 + ^glaF^]Fi 

An 
+ 45 (1 - —)F^}/K2 - (a - a')/(272), 

72 

(38.3) 

(flta) = 2[(a - a'fF2 - 4,gfa{a - a')^3 

+ 2gfF,]/K2, (38.4) 

1 



(-Raa^a) 



(k + 71)^272 



-2(a - a'YxoF2 



+ 8gla{a ~ a')xoF3 - Agl{-/^ + 4a(a - a'))F5 

-4g2«7|((«; + 27i)2-4aa')i^3 
-{a- a'fK2l2] 
2.91 



(38.5) 



{R3a)s = --^[-4(k7i + 272a)aa' 
i<i 

+ 272 (k^ + 47I + 4K7i)a 

+ k{kji + 2jI){2ji + k)], 

2 



{R3a')s = 



{[{a - a)xiF2 - Ag^axiFz 



2glx2F5\lK2 + X3l^}, 



(38.6) 



(38.7) 



Ki = {[Aaa' - k(27i + K)f - 167^aa'}7i, 
K2 = {256(aa')'* - 64(4^2 + 6K71 + b-il){aa'f 
+ 16[2(k2 + 2^71 + 27|)(2k2 + 2^71 + 7I) 
+ n\{2^,^nf + {n + ^,f)]{aa'f 
- 4k2[2(k2 + 2k7i + 27^) (71 + nf 
+ (2^2 + 2^71 + 7|)(k + 27i)2]aa' 
+ k4(k + 7i)2(k + 27i)2}7i, 
K3 = [Aaa' — k(7i + k)] — 4j2CtOi' , 



with 



2^0 = 72 + '*('* + 71) ~ 4qq;', 
xi = Xo(a + a')/72 + '«('« + 7i), 
X2 = -72 + 4a(Q: + a')/72: 
X3 = (a^ - a'^)/72 

- 8g?K72[(K + 27i)2 - 4aa']K3/K2, 

" = *2f]' 

" =* — 2n — ■ 



where 



Fi = — 4aa'7i - 272^0' + (271 + k)(k7i + 272), 
F2 = 16(aa')^7i 

- 4[37i(a. + 7i)2 + 27f + (7I - 72)(^ + 57i)](aa')' 
+ {7iK^(27i + nf + 2(k + 7i)(k7i + -fi)[(K + -fif 
+ 27|-7?]}a«' 

- «;2(k + 7i)(k7i + 7l)(27i + k)V4, 
F3 = 32(aa')^7i 

- 8[(k7i + 27|)(27i + k) + {k + ji){kji + ji)]aa' 
+ (271 + k)(7ik + 27^)(2k^ + 2^71 + 7^) 

+ 72(271 + «)'(7i + 2k) - 472(71 + K){-fl - 72), 
Fi = [k(k + 71) - 4aQ;']/72 + a + a\ 
F5 = -64(aa')^7i 

+ 16[(3k7i + 2-iI)(k + 71) + 71(372' + K-ii)]{aa'f 

- 4[(k7i + 27|)(27i + k)(2k7i + 2k' 
+ 7ik2(k + 7i)']aQ;' 
+ «;2(27i + «;)(«;7i + 27|)(«; + 71)2, 



-7I) 



APPENDIX II 



In this appendix, we present the Laplace transform of the two- 
time correlation function of the cavity field operators appear- 
ing in the expression for the incoherent part of the spectrum. 
We have applied the quantum regression theorem to Eqs. (|20l i 
and (ISTT i and calculated the Laplace transform with the ini- 
tial condition given by the stationary solutions listed in the 
Appendix I for the case of the cavity frequency tuned on reso- 
nance with the central frequency of the dressed-atom system. 



(a'(p),a)s = ^, , ,^, , , (39) 



2(p + 7i)Z?(p) 



where 



D{j>) = aa'^l - [{p + |)(p + I + 71) - aaf, (40) 
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and 

mo = [2{gi{R3a)s + 71 (a^a).) (71 + «/2) + {{a - a'){R3a^), 
+ {a + Q;')(i?3a^a)s)7i][-K(7i + k/2)/2 + aa'] 
+ 27i[k(7i + k/2) - 2aa'](7i + K){a)s{a^)s 

+ 72(71 + K/2)[7i((a + a')(a^a)s + (" - a')(a^)s) 
+ 2gia{R3a)s] + 2jij2aa' {Rsa^ a) s 

+ 251(71 (72 - a)(7i + k) + 72'«(7i - '«^/4 + 7i'«/2 
+ 72a + aa')/(27i)](a)s, 

TOi = [2gi(i?3a)^ + (471 + K){a^a), + {a - a')(^3a^)s 

+ (a + a'){R3a'^a}s][-Ki-fi + k/2)/2 + aa'] 

- (71 + K)[2(gi{R3a)s + -fi{a''a)s){ji + k/2) 
+ 71 (a - a'){R3a^)s +71 (a + a' ) {R^a^ a) s] 

+ 2[7i((7i + nf + k(7i + k/2) - 2aa') + (^(71 + k/2) 

'^2aa'){-ii + n)]{a)s{a))s 

+ [(271 + K/2)((a + a'){a)a)s + {a - a){a^)s) 

+ 2gia{R3a)s + 2aa' {R^a^ a) s]j2 
+ 25i[(72 - a)(27i + k) + 72K/2](a),, 
m2 = -[{a- a'){R3a'^)s + (a + a')(i?3a^a),](K + 271) 
-gi(47i+3K)(i?3a), 
+ [-71(471 + 3k) - (3k/2 + 7i)(k + 271) 
+ 72(0: + a') + 2Q!a'](a^a)s + 72(0; - a'){a^)s 
+ 2[(7i + nf + 271(71 + k) + k(7i + k/2) 

- 2aa]{a)s{a^)s + 2.91(72 - ")(«)«, 

"^3 = -(a - a')(^3a^)s - (a + Q:')(i?3a'l'a)s - 2gi(i?3a)s 

- 3(271 + n){a)a)s + 2(371 + 2K)(a),(at),, 
TO4 = 2((a)s(at)5 - (ata)s). 



